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We present an effective theory for the Bloch functions of a two-dimensional square lattice near 
a quadratic degeneracy point. The degeneracy is protected by the symmetries of the crystal, and 
breaking these symmetries can either open a bandgap or split the degeneracy into a pair of linear 
degeneracies that are continuable to Dirac points. A degeneracy of this type occurs between the 
second and third TM bands of a photonic crystal formed by a square lattice of dielectric rods. We 
show that the theory agrees with numerically computed photonic bandstructures, and yields the 
correct Chern numbers induced by parity breaking. 



In a two-dimensional crystal with a square-lattice 
{Civ) symmetry, it may happen that a pair of bands are 
degenerate at a point of high symmetry, such as the cen- 
ter (r) or corner (M) of the Brillouin zone. An example 
of this occurs in a photonic crystal formed by a square 
lattice of dielectric rods. As shown in Fig.flja), the sec- 
ond and third TM bands exhibit a quadratic degeneracy 
at the M-point. The existence of this degeneracy is inde- 
pendent of details such as the permittivity and radius of 
the rods, as long as the C± v symmetry is preserved. This 
degeneracy is of particular interest since Wang et. al^ 
have recently shown that it can be lifted by gyromagnctic 
effects, which break parity and time-reversal symmetry, 
and that the bandgap opened in this way is populated 
by a one-way edge mode analogous to chiral electronic 
edge states in the quantum Hall effect. Electromagnetic 
one-way edge modes were first predicted by Haldane and 
Raghu^, who argued that they typically occur in sys- 
tems possessing "Dirac points" , meaning that the modes 
near each degeneracy point can be described by an effec- 
tive Dirac Hamiltonian. Although Dirac points have been 
extensively analyzed in the condensed-matter literature^, 
there has been, to the best of our knowledge, no analo- 
gous study of these quadratic degeneracies. Since the 
system proposed by Wang et. al. currently appears to be 
the most promising for realizing electromagnetic one-way 
edge modes, due to its large relative bandgapi, there is a 
present need for an effective theory of such degeneracies. 

In this paper, we present an effective theory that de- 
scribes the bands near a degeneracy point, based on the 
symmetry properties of fc-space around that point. We 
show that the quadratic degeneracy in the C± v crystal 
can be regarded as a pair of linear degeneracies, ana- 
lytically continuable to Dirac points, that are "pinned" 
to the same fc-space point by the crystal symmetry. 
The quadratic degeneracy is robust against perturba- 
tions that preserve this symmetry. It can be lifted by 
parity and time-reversal symmetry breaking (which we 
will henceforth simply refer to as parity breaking.) In 
that case, the two bands acquire Chern numbers^^ of 
±1, in agreement with the numerical result of Wang 
et. alX Breaking the 90° rotational symmetry "unpins" 
the quadratic degeneracy point, which splits apart into 
two distinct linear degeneracies. The theory applies to 
any two-dimensional Bloch system, whether electronic or 
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FIG. 1: (Color online) TM bandstructure of a two- 
dimensional photonic crystal formed by a square lattice of 
dielectric rods, (a) Fully symmetric crystal with rod radius 
0.15a (where a is the lattice constant) and e = 20. The rods 
are embedded in air, and /j, — 1 everywhere. The quadratic 
degeneracy between the second and third TM bands is indi- 
cated with an arrow. Inset: the crystal structure in real space, 
(b) Parity-broken crystal, with off-diagonal permeability com- 
ponent fi xy — O.li. Inset: Band structure near the lifted de- 
generacy, with dots showing numerical solutions of the exact 
Maxwell equations^. The solid lines show the analytic ap- 
proximation of Eq. [5j where Ao, P, and 7 are calculated from 
the symmetric system and Q2 is calculated separately based 
on the proportionality constant with \fi X y\ (see Fig. [4]) 



photonic, with C^ v symmetry and a quadratic degeneracy 
point. In particular, wc show that it accurately describes 
the aforementioned photonic crystal of dielectric rods for 
a wide range of dielectric contrasts and rod radii. 

Let us consider a crystal in which two bands are degen- 
erate at a point k = km, with unbroken C^ v symmetry. 
We choose a pair of independent Bloch functions at /cm, 
denoted by u x M {r) and u 2 M {r). The Bloch functions at 
neighboring values of k can be written as 

u\{r) = Cu(fc)n M (r) + c 12 {k)u 2 M {r) 
u\{r) = c 21 {k)u x M {r) + c 22 {k)u 2 M {r). 

The mixing elements c nm {k) can be related to the mode 
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frequencies, u) n (k), through an eigenvalue equation 



H(k) 



Cnl (k) 



Cnl (A) 
Cn2{k) 
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where the "effective Hamiltonian" H (k) is a 2 x 2 matrix 
whose eigenvalues are, by definition, X n (k) = uj n (k) — luq. 
We will not be concerned with the value of the "zero- 
point" frequency lu . Now, suppose we alter the system 
that we used for defining u l M {r) and u 2 M {r), such as by 
breaking some of its symmetries. If the perturbation is 
sufficiently weak, the Bloch functions of the altered sys- 
tem can still be described by |T]) for k ~ fcjwj with some 
new choice of c nm (k) and hence of H(k). 

In this system, we will be interested in three different 
symmetry-breaking operations. Firstly, we could "shear" 
the lattice by rotating the basis vectors as follows: 



a\ = a (cos 9, sin ( 
ti2 = a (sin 9, cos ( 



(3) 



where a denotes the lattice constant. This breaks the 
symmetry under C4 rotations and reflections about the 
x and y axes. Secondly, we could distort the rods by 
stretching them along the x or y axes — or, alternatively, 
stretching the lattice vectors and rescaling k x and/or k y ; 
this breaks the symmetry under rotations and reflections 
about y = ±x. Thirdly, we could break parity, which 
can be accomplished in a photonic crystal, e.g., using a 
magneto-optic effect that adds an imaginary off-diagonal 
term \i xy = irj to the permeability tensor of the rods^. 

The goal is to find an effective Hamiltonian that de- 
scribes bands such as those in Fig. [TJ including the results 
of the above symmetry-breaking operations. We claim 
that the desired Hamiltonian has the following form: 



H= A r 



E 
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FIG. 2: (Color online) Photonic bandstructures for time- 
reversible lattices of dielectric rods with radius ro and per- 
mittivity e, plotted against the fc-space displacement from 
the corner of the Brillouin zone, it = k — (n/a, 7r/o), along the 
line k — (k x , —Kx)- Dots show numerical data^; the solid lines 
show the analytic results after fitting the free parameters in 
Eq. to the numerical data. Note that the bandstructure is 
independent of the parameter /3 along this given line. In (a) 
and (b), we take ro = 0.25a and e = 16.26, with fitted value 
7 ~ 0. For (a), 6 = (square lattice) and «3 = 0; for (b), 
9 = 1.2xl0" 4 radians with fitted value 013 = 3.7x 10" 5 (ir/af. 
In (c), (d), and (e), we take ro = 0.2a and e = 9.92, with fitted 
value 7 ~ -0.5. In (c), = and a 3 = 0; for (d), 6 = 10~ 4 
radians with fitted value q?3 = 1.8 x 10~ 5 (-7r/a) 2 . In (e), a 3D 
plot of the bandstructure of system (d) is shown. 



where k = k — Um is the fc-space displacement from the 
degeneracy point and Oi are the usual Pauli matrices. 
The phcnomcnological parameter Ao determines the fre- 
quency scale, (3 and 7 control the relative curvatures of 
bands along different directions, a\ controls the relative 
lengths of the two lattice vectors, 011 is proportional to 
the parity-breaking permeability component 77, and 03 is 
proportional to the shear angle 9 defined in This 
Hamiltonian is valid in the neighborhood of k = 0, and 
we have omitted 0(n 4 ) terms which have negligible ef- 
fects on the band properties in the regime of interest. 
Furthermore, we assume that the symmetry-breaking is 
weak (e.g. 9 -C 1), and thus retain only symmetry- 
breaking terms that are zeroth-order in n. 

Before discussing the validity of the ansatz let us 
determine its band structure and then show that it is con- 
sistent with the symmetries of the system. First, consider 



a.\ = «3 = 0, for which the eigenvalues are 
A±(/S)/A = l\A 2 ± v /|k| 4 + (/3 2 -1)(^-^)2 + q 



2 

2 ■ 

(5) 

Suppose we assume (3=1 (setting /3 ^ 1 simply dis- 
torts the bandstructure along the k x = ±n y directions). 
For oli = 0, (O then reduces to a pair of quadratic bands 
A±/Ao = (7±1)|k| 2 , which meet at k = 0. The parameter 
7 controls the relative curvatures of the two bands. For 
instance, when 7 = 0, the bands have equal and opposite 
curvatures, as shown in Fig. Ela); when —1 < 7 < 0, 
the two bands curve in opposite directions but the upper 
band is flatter, as in Fig. [UJc). Setting a 2 7^ lifts the 
degeneracy and opens a bandgap AA = 2Aoai2- The two 
bands will curve in the same direction at k = 0, as ob- 
served in Fig. 2 of Wang et. alX and in Fig. HJb) of the 
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current paper. 

Next, let us consider a 3 7^ 0, keeping a\ = 0. For 
Q2 = 0, the quadratic degeneracy splits into two dis- 
tinct degeneracy points, at k± = ±( y/ a 3 /2, — y/a^/2 ) 
for a 3 > and K± = ±( ^/— a 3 /2, a 3 /2 ) for a 3 < 0. 
We will henceforth assume that 013 > 0; the following dis- 
cussion can be easily adapted to the a 3 < case. Let us 
expand the Hamiltonian around k± , using the variables 

Si = o y Kx + K v> 

I (6) 

which arc simply fc-space displacements from k± , rotated 
by 45°. To first order in q\ and q-2, 

H±(q)/X Q ~ 7(a 3 =F V^Sa^ g 2 ) 

± a/803 (PqiCTl + 920" 3 ) + a 2 cr 2 . 

When 7 = and = 1, this reduces to a two-dimensional 
Dirac Hamiltonian near each each degeneracy point (or 
"Dirac point"). Furthermore, a 2 plays the role of a mass 
term, opening a bandgap AA = 2Ao« 2 . Setting 7 7^ 
distorts the Dirac Hamiltonian and its eigenvalue spec- 
trum; for example, the Dirac cones in the a 2 — limit 
become tilted in fc-space, as shown in Fig. Eld). Along 
the line n y = ±k x , the splitting of the degeneracy point 
can be thought of as a vertical relative displacement of 
the two parabolic bands (note, however, that the bands 
meet only at isolated points in the full /c-space.) The 
splitting is accompanied by a change in the density of 
states from a discontinuity to a linear "dip" centered at 
the frequency of the band degeneracy. When a 2 ^ 0, 
the density of (TM) states is discontinuous at the band 
edges, dropping to zero inside the band gap. 

The situation is very similar for ol\ 7^ 0. When 
Q-2 = «3 = 0, the degeneracy splits into two, but along 
the line K y = (if ai > 0) or n x = (if ai < 0), in- 
stead of K x = ±K y . When both ot\ and 03 are nonzero, 
the degeneracies are located at an intermediate location, 
k± = ±((2j + a 3 ) 17,4 (cos (j>, sin 0), where tan</> = a 3 /ai, 
and expanding around each point yields a Dirac-likc 
Hamiltonian analogous to ([7]). 

When a 2 7^ 0, the bands are non-degenerate, and their 
Chern numbers^ can be calculated. The details of this 
calculation are given in the Appendix, and the result is 
that the upper and lower bands possess Chern numbers 
— sgn(a 2 ) and sgn(a 2 ) respectively, regardless of the val- 
ues of ax, 03, /3, and 7. This implies the existence of a 
single family of one-way edge modest, and agrees exactly 
with the numerical results of Wang et. alX. Although 
the effective Hamiltonian (U|) is only valid near k = 0, 
it yields the same Chern number as the actual band- 
structure because only the region near the broken degen- 
eracy point provides a non-vanishing "Berry flux" con- 
tribution to the Chern number— Furthermore, while 
our theory only describes weak symmetry-breaking, the 
Chern number is a topological quantity and cannot be 



altered by non-perturbative distortions (as long as the 
bands remain non-degenerate^), which is why it remains 
unchanged even in the strong parity-breaking regime ex- 
plored by Wang et. al. When a\ and/or 0:3 are non-zero, 
the two linear degeneracy points each contribute ±1/2 
to the Chern number, in accordance with previous anal- 
yses of the Dirac Hamiltonian^. When a\ = a 3 = 0, the 
Berry connection winds twice as fast around the point 
K = 0, which provides the entire contribution of ±1. The 
dependence of the Chern number on the sign of a 2 con- 
firms that a 2 controls parity breaking, since the Chern 
number can be shown to vanish identically when parity 
is unbroken. 

The fully symmetric Hamiltonian Hq = H\ ai —o must 
transform under any operation g G C^ v as 

D(g)H (it)D-\g) = H (gH). (8) 

Regardless of the values of (3, 7, and Ao, this holds if 
D(g) falls under E, the only two-dimensional irreducible 
representation of C^ v . Thus, ±90° rotations can be rep- 
resented by ^i<J2, reflections about the k x (k v ) axes by 
<7i (—01), and reflections about K y = ±k x by ±er 3 . 

By studying how H(iZ) transforms under E, we can 
show that the quadratic degeneracy is protected by the 
crystal symmetry. Any zeroth-order term proportional 
to the identity matrix, when added to Hq, simply shifts 
the eigenvalues without opening a gap. Adding a zeroth- 
order term proportional to a\ (i.e. a± 7^ 0) breaks C± v 
since, under the representation E, a — > —a for 90° rota- 
tions and reflections across k x = ±K y . Note that a — > a 
for reflections across the k x and n y axes, in agreement 
with our claim that a\ 7^ corresponds to stretching the 
lattice. Similarly, setting a 2 7^ preserves the rotational 
symmetries but breaks the reflection symmetries (parity) . 
Finally, setting a 3 7^ preserves the reflection symme- 
try across k x = ±n y but breaks the symmetry under 90° 
rotations and reflections across k x = and K y = 0. 

Furthermore, the Hamiltonian cannot include terms 
that are first-order in k if the C± v symmetry is unbroken 
or only partially broken. Such terms have the general 
form AH = J2i=i Y^=i K i c ij a ji an d we can show that 
Cij = for all i, j as long as the system is symmetric un- 
der either rotations, reflections about the k x and thy clXGS, 
or reflections about k. x — ±n y . At least one of these sym- 
metries is preserved by each of the three "elementary" 
distortions discussed above. This situation may be con- 
trasted with a triangular or honeycomb lattice, for which 
there is a C 3lI symmetry around each corner of the hexag- 
onal Brillouin zone. There, one can write down an O(k) 
Hamiltonian which transforms under a two-dimensional 
irreducible representation of C 3lI : this is just the Dirac 
Hamiltonian 

H'(k) = A (kztri + k y cr 3 ) . (9) 

In this case, a zeroth-order "mass" term proportional to 
02 controls parity breaking^. 

We checked the validity of the ansatz Q and com- 
puted the parameters o^, /3, and 7 using the MPB^ and 
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FIG. 3: (Color online) Values of 7 and j3 obtained from least- 
squares fits of Eq. ([5]l to bandstructures computed numeri- 



cally along the lines k y = n/a and k x 



The photonic 



crystal consists of a square lattice of rods with radius ro and 
permittivity e, embedded in air, with /j, = 1 everywhere. 



COMSOL2 computer programs, which solve the Maxwell 
equations without any approximations apart from the 
discretization of the simulation cell. We first determined 
(3 and 7 for ot\ = U2 — &3 = by fitting the computed 
bandstructures of the fully-symmetric crystal to ©. As 
shown in Fig. [3], we obtain values for j3 and 7 that depend 
on the details of the crytal, including the rod radius and 
permittivity. 

Upon perturbing the simulated crystals by changing 
the lattice constant slightly along (say) the x direction, 
8a x < 0, the degeneracy splits along the line n y = in 
the computed bandstructures. From the location of the 
linear degeneracies, we can obtain a%, which turns out 
to be proportional to 6a x . Similarly, setting 9 7^ in- 
duces linear degeneracies along the line k x = K y , and 
a3 is found to be proportional to 9. By computing the 
bands for off-diagonal permeability component fj, xy = ir\ 
in the rods, and fitting these to ((SJ), we find that 0L2 is 
proportional to r\. These results are shown in Fig. 3] 
Like (3 and 7, the proportionality factors a\/8a x , 0*2/11, 
and CV3 /9 depend on the details of the crystal. The fact 
that the linear degeneracies induced by ^ are not 
determined solely by the lattice geometry stands in con- 
trast to Dirac points in previously-studied triangular and 
hexagonal lattices, which are pinned to the corners of the 
Brillouin zone. 

In conclusion, we have presented an effective theory 
for a pair of bands in the vicinity of a band degeneracy, 
within a two-dimensional C^-symmetric crystal. The 
high degree of symmetry near the degeneracy point deter- 
mines the form of the 2x2 effective Hamiltonian matrix, 
including the lowest-order symmetry- breaking terms. In 
particular, we showed that quadratic nature of the dc- 
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FIG. 4: (Color online) Plots of ai, 02, and 03 for two dif- 
ferent photonic crystals. In the left column, the dielectric 
rods have radius ro = 0.1a (where a is the lattice constant) 
and permittivity e = 15. In the right column, ro — 0.3a 
and e = 20. We obtain a\ from the location of the de- 
generacies k = ±(0, y/ai/P) for lattice stretching parame- 
ter 8a x /a. We obtain 0:2 by performing a nonlinear least- 
squares fit of Eq. ([5J to the computed T-broken bandstruc- 



tures along the line k x = 



using the values of 7 found 



in Fig. [3] (which were obtained using the symmetric lattice); 
it is plotted against the imaginary off-diagonal permeability 
component, fj, xy = irj. We obtain 03 from the location of the 
degeneracies k = ±(^/ce:}/2, —y/az/2"), for lattice distortion 
angle 9. 



generacy is protected by the C^ v symmetry. We have 
also shown that the theory accurately describes the sec- 
ond and third TM bands of square-lattice photonic crys- 
tals near their M-point degeneracy, including the correct 
Chcrn numbers when parity is broker^. On the other 
hand, the theory should be applicable to any electronic 
or photonic system possessing C^ v symmetry and a two- 
fold band degeneracy; it is possible, for instance, to con- 
struct tight-binding electronic models that reduce to our 
effective Hamiltonian (jj) near a degeneracy point. The 
advantage of the present method is that it relics only 
on symmetry principles, and can therefore be applied to 
systems, such as photonic crystals, where other meth- 
ods such as tight-binding are not applicable. For lattice 
symmetries other than C± v , effective Hamiltonians can 
be constructed by finding an appropriate representation 
of the symmetry group. 
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APPENDIX: CALCULATING THE CHERN 
NUMBER 



The integral can be performed via the substitution 
sin 20 = tanhu, and we obtain 



In this Appendix, we describe the calculation of the 
Berry connection and Chern number— for the bands asso- 
ciated with our effective Hamiltonian. We will consider 
the lower band the calculation for the upper 

band proceeds analogously. 

First, consider ot\ = 0:3 = 0. We note that the eigen- 
vectors of the effective Hamiltonian (|4]) do not depend 
on 7 since that parameter multiplies the identity matrix. 
For simplicity, we set (3 = 1 . The eigenvector correspond- 
ing to the lower band is 



IV--(k)} = 



1 



^2[n A + al + K 2 y^FTa^ sin 20] 
— n 2 cos 20 + i/3 



k 2 sin 20 + ^k 4 + i 



(A.l) 



regardless of the values of 7. Here, (n, 0) is the cylindrical 
coordinate representation of k. The Berry connection is 



A-(k) = ^-(k)\V k ^-(k)) 

ict2K ^cos 20 k — sin 2(f) 
k 4 + /?V k4 + «2 sin 20 



(A.2) 



To obtain the Chern number, we integrate the Berry con- 
nection around a loop k = kq: 



C = — <b dk ■ A~(k) 
2?" J k=ko 



2a, 



k,q sin 20 1 



(A.3) 



K Q + Q 2 + K Q \/ K + a 2 Sm 20 



c- 



sgn(a 2 ) 



a 2 



(A.4) 



sgn(a2) for |a 2 | <€. Kq- 



As discussed in the text, the above result remains un- 
changed even when we enter the non-perturbative regime, 
even though our effective theory is only valid for small 
values of k and «j. 

When a.\ and/or are non-zero, the band maximum 
at k = splits into two distinct maxima, and expanding 
around each maximum yields a Dirac-likc Hamiltonian. 
For instance, when ot\ = and 03 7^ the maxima occur 
at k± = ±( \J as /2, — sJa~sj2 ), and the Hamiltonian near 
each of these points is given by ([7]). In terms of the 
variables q\ and <?2 defined in ([5]), the Berry connection 
for the lower band is 



t* ib 



cos< 



sine 



■ — , (A.5) 
q 2 + b 2 ± q\/q 2 + b 2 sin 
where 6 = ^2/ y/ots , ± refers to which maximum we are 
expanding around, and (q, 0) is the cylindrical coordinate 
representation of q. This Berry connection has the same 
form as (|A.2j) . but winds half as quickly around each 
maximum point as (|A.2[) docs around k = 0. Each max- 
imum thus contributes sgn(a 2 )/2 to the Chern number 
of the lower band. 
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